(N 

o 

(N 
Of 

On! 



a- 

1). 

fi: 

(N : 

> 
(N 

in 

ON' 

O! 
(N 



Stability of Leptonic Self-complementarity 
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We investigate a stability ol leptonic self-complementarity such that sum of three mixing angles in 
lepton sector is 90 degrees. Renormalization group equations in a context of minimal supersymmetric 
standard model for the self-complementarity are analyzed. It is seen that one of Majorana phases 
plays an important role for the stability of self-complementarity. We find some stable solutions 
against quantum corrections at a low energy. An effective neutrino mass for neutrino-less double 
beta decay is also evaluated by the use of neutrino parameters giving rise to the stable solutions. 

PACS numbers: 14.60.Pq, 12.60.Jv 



Neutrino oscillation experiments established that there 
are two large mixing angles (6*12 and $23) of Pontecorvo- 
Maki-Nakagawa-Sakata (PMNS) matrix in lepton sec- 
tor. Then a non-vanishing #13 in the PMNS has been 
reported by recent long baseline and reactor neutrino 
experiments [l|. These results can be interpreted by 
three flavor mixing of neutrinos. Regarding with neu- 
trino masses rrii (i = 1, 2, 3), the neutrino oscillation ex- 
periments determine only two mass squared differences, 
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|m 2 | - an d lAm^l 
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Therefore, two types of neutrino mass hierarchy are al- 
lowed, i.e. normal hierarchy (NH) mi < 1112 < and 
inverted hierarchy (IH) 7713 < mi < m-j.. Further, neu- 
trino experiments have not determined whether the neu- 
trinos are Dirac or Majorana particles. Clearly, the na- 
ture of neutrinos would be a key to find physics beyond 
the standard model (SM). 

In theoretical side of neutrino physics, various ap- 
proaches have been discussed in order to investigate hid- 
den flavor structure behind the SM, e.g. introductions 
of flavor symmetry, mass (matrix) texture analyses, and 
searches for exotic relations among flavor mixing an- 
gles etc.. In this work, we focus on a leptonic self- 



complementarity (see also @ for related discussions) 
as 
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(1) 



We start with effective Yukawa interaction and Wein- 
berg operator at a low energy scale such as electroweak 
(EW) scale Aew in a context of minimal supersym- 
metric standard model (MSSM), Ly = —y e LzHdeji + 
^(H u L L )(H u L L ) + h.c, where L L are left-handed lepton 
doublets, eji are right-handed charged leptons, H U {H^) 
is up(down)-typc Higgs, y e is Yukawa matrix of charged 
leptons, and k(H u Lt){H u Ll) is the Weinberg operator, 
which can be effectively induced by integrating out a 
heavy particlc(s). One of examples to obtain this op- 
erator is seesaw mechanism. Typical scale of the see- 
saw mechanism is O(10 14 ) GeV. Therefore, note that the 
effective coupling k is having mass dimension —1 and 
k _1 ~ O(10 14 ) GeV. Such a heavy mass scale can real- 
ize tiny active neutrino mass scales through the seesaw 
mechanism. In this work, we utilize an useful parameter- 
ization for the PMNS matrix as 
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where = sin0y, c = cos0jj (i,j = 1,2,3; i < j), S is 
a Dirac phase, and D p is a diagonal phase matrix includ- 
ing two Majorana phases, p and a. An neutrino mass 
matrix M v can be diagonalized as Vj M U V* = M^ iag = 
Diagjmi, m2, 7713} with rrii = Kiv\ where v u is vacuum 
expectation value of up-type Higgs. 

Next, we consider renormalization group equations 
(RGEs) in the MSSM. The RGE of k is given by 
16-7T 2 (dn j dt) — an + [(y e yl)n + n(y e y\) T ] with a = 
6[— ffi/5 — I? 2 + ^{UuVu)] where gi are gauge coupling 



constants, t is an arbitrary renormalization scale as t = 
ln(/z/Ab and A is a high energy scale such as the seesaw 
scale [a, @] • One can also obtain RGEs of 0y in a diagonal 
basis of y e as ddij/dt = ^(012,023,013, Ki,S, p,a,y T ;t) 
where right-hand side (RHS) of this equation is given 
in Q (see also Q for other discussions of mixing an- 
gles under the RGEs). Now we turn to the self- 
complementarity relation (JTJ) and investigate the follow- 
ing equation, 
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fa^2° l 3 = X] F v(^2, 6,23 ' 6,13 :Ki,6,p,<r,y T ;t) 

ij ij 

= F(6 12 , 6*23, 6*13, «tj *j P, cr, j/ T ; i), (3) 

where ij is summed over 12, 23, and 13. The function F 
is described by 3 mixing angles, 3 effective couplings for 
the light neutrino masses (or equivalently light neutrino 
masses m,), 3 CP-phases, a Yukawa coupling of t, and 
renormalization scale. Then once we impose (fTJ) on ([3]) at 
an energy scale to, one of mixing angles in F is removed 
as e.g. F(6i2, 023, Ki,S, p, a, y T ;to). We now focus on an 
equation, 

F(8i2,023,Ki,S,p,a,y T ;t ) = 0. (4) 

This equation means that once the equation is satis- 
fied at an energy scale to, the self-complementarity is 
also satisfied at all other energy scales t, i.e. the self- 
complementarity is stable against quantum corrections, 
if running effects of parameters except for mixing angles 
are tiny. In fact, we can find consistent solutions of (|4]) 
with experiments for both NH and IH cases. According 
to the latest experimental data of neutrino oscillation [8[ 

31.3° < flu < 37.5°, (5) 
38.6° < 6*23 < 53.1°, (6) 
7.0°(7.3°) < flu < 10.9°(11.1°), (7) 

at 3a level for the NH(IH), the (QJ can be satisfied. 

Mass spectra of neutrinos at a low energy are 
defined by mi = yf m\ — | Am§ x | and = 
\Jm\ — |Am|i| + Am 2 ! with best fit values Am|i = 
7.62 x 10~ 5 eV 2 and Am^ = 2.53 x 1CT 3 eV 2 for 
the NH, and mi = \J m\ — Amji and TO 3 = 
y/m\ - \Am\ Y \ - Am^ with AmJi = 2.40 x 10~ 3 eV 2 
for the IH. Therefore, the largest neutrino mass 7713(7712), 
of NH(IH) case is a free parameter in our analyses. 
We analyze in range of y/\ Am| x < 7773 < 0.2 eV 
(y/lAm^l + Ajti 2 ! < m 2 < 0.2 eV). The case of 
7713(777,2) = 0.2 eV corresponds to a degenerate mass spec- 
trum. In such case, 777.3(771.2) is bounded by a cosmological 
constraint on the sum of neutrino mass as m^ < 0.6 
eV P, and thus as m 3 (m 2 ) < J2 m i/3 — 0.2 eV. There- 
fore, this must be implied as upper bound on the largest 
neutrino mass. Here note that since input values of neu- 
trino parameters at a low energy are used, the solutions 
of are corresponding to ones at to ~ 1ii(Aew/A). 
The y T has been also approximated at a low energy as 
j/ t (Aew) = 10~ 2 in our analyses. 

There are 7 parameters (3 mixing angles, 3 CP-phases, 
and 1 neutrino mass, 7773 or 777,2) and 2 imposed equations 
(([1]) and ((4])). Therefore, number of free parameters is 
5. Since it is however intricate to deal with all 5 pa- 
rameters as completely free ones, we numerically analyze 
at some fixed neutrino masses as examples. According 



(a) m 3 = y^Amfj ( b ) m :i = y^jAmfJ + Am| 




FIG. 1: Examples of solutions of ([4]) 



to our analyses, a CP-phase is important to give solu- 
tions of (J4j) with ([l]); we numerically found that there is 
no solution to satisfy (0| in cases of (5 = p = a = 0), 
(5^0, p = a = 0), (p^ 0, 6 = a = 0) for both NH and 
IH but (a ^ 0, 6 = p = 0) can give solutions of ((4]). In 
the following, we focus on the case of (er ^ 0, 5 = p = 0) 
in detail. Now we have 2 free parameters (a and one of 
mixing angles) in order to look for solution, i.e. once we 
fix a and one of mixing angles, all values of our parame- 
ters are uniquely determined as we will explain below. 

We have scanned over < \a\ < tt. Some results of 
numerical analyses are shown in FIG. [1] as examples. 
In the figures, vertical and horizontal axes are 023 and 
012, respectively, and shaded region mean that the self- 
complementarity (fTJ) is correlatively satisfied within 3ct 
ranges of mixing angles ([5]l-(|7]l. Lower(upper) slanting 
and right(lcft) sides are bounded by maximal (minimal) 
013 and 0i2 at 3a level, respectively. Note that an allowed 
range of 023 becomes narrow as 41.6° < 023 < 51.7° com- 
pared to (JHJ) due to (fTJ). Both two lines in the figures 
are solutions of (U]), which are contours of a. FIG. [1] (a) 
is case of 7773 = ^/|Am§i| for the NH where mi = 0, 
and there are 2 remaining parameters, a and one of 
mixing angles. In the case, 2 lines corresponding to 
er| ~ 126.0° and 132.0° arc grazing shaded region at 
upper left described by A and lower right B points, re- 
spectively (there are 2 (different sign) solutions of a in all 
cases, i.e. solutions are symmetric for reflection respect 
with 7r). Therefore, values of all 7 parameters in our anal- 
yses are uniquely determined at e.g. point A or B. Fur- 
ther, the points A and B in this case determine maximal 
and minimal values of parameters for the stability of self- 
complementarity. Therefore, the self-complementarity 
is stable in shaded regions of 126.0° < \a\ < 132.0°, 
31-3° < 0i2 < 37.5°, 7.0° < 0i3 < 10.9°, and 
Min[0 23 ] < 023 < Max[0 23 ] where Max(Min) [0y] is maxi- 
mal(minimal) value of Otj. Max(Min)[023] is evaluated by 
Max(Min)[0 23 ] = 90° - Min(Max)[0i 2 ] - Min(Max)[0i 3 ] 
due to |T]). Note that value of one of mixing angle (e.g. 
023) is not independently taken because of ([T]). In the 
case, we obtain Max(Min)[0 23 ] ~ 51.7°(41.6°). The re- 
sults are summarized in TAB. HI 
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NH 


IH 


777-3 




V\Am 2 31 \ + Am^ 


0.2 eV 


eV 


y/ml- lAmfjl - Am^j 


777,2 


V Am 21 


1.23 x 10~ 2 eV 


0.194 eV 


v/[Am^l + Am^ 


5.73 x 10" 2 eV 


777 1 


eV 


8.73 x 10~ 3 eV 


0.194 eV 


V\ Am ii\ 


V m 2 ~ Am 21 


JVLm[0i2(sm P12JJ 


31.3° (0.27) 


31.3° (0.27) 


Max|#i 2 (sm 6^12 JJ 


37.5° (0.37) 


37.5° (0.37) 


31.5° (0.27) 


Min[0 2 3(sm #23)] 


41.6° (0.44) 


41.4° (0.44) 


47.4° (0.54) 


iviax[t/23 (sm P23 J J 


51.7° (0.62) 


50.7° (0.60) 


51.7° (0.62) 


51.4° (0.61) 


48.0° (0.55) 


Min 173, „f sin' 4 ft, o)l 


7.0° (0.015) 


7.3° (0.016) 


7.0°(0.015) 


7.3°(0.016) 


10.7° (0.034) 


lVl<AX|yi3 ^blll C13 


10.9° (0.036) 


11.1° (0.037) 


Mill [cr] 


126.0° 


124.9° 


102° 


78.8° 


87.6° 


Max [cr] 


132.0° 




105° 


80.4° 


88.8° 


Min[(m) ee ] [meV] 


2.36 


3.61 


73.5 


14.9 


25.7 


Max[(m) ee ] [mcV] 


3.09 


6.83 


97.4 


23.1 


26.1 



TABLE I: Examples of solutions, and minimal and maximal values of neutrino parameters in the corresponding regions 



FIG. [T] (b) shows a case of m 3 = y/\ Am§ x | + Am^. 
In the case, contour of \a\ — 7T, which is the turn of con- 
tours of solutions, is across the shaded region and the 
value of \a\ — tt is a maximal value of |oj. Therefore, 
solutions can be obtained in range of 124.9° < |ct| < 
tt. We find that points A-C determine Max(Min)[0y] 
as 31.3° < 6»i2 < 37.5°, 7.3° < 6 13 < 10.9°, and 
Max(Min)[6> 23 ] ^ 50.7°(41.6°). Derivations of Max[6» 23 ] 
and Min[#i 3 ] are more complicated than the previ- 
ous case. Max [6*23] in this case is not determined by 
Max[(9 23 ] = 90° - Min[(9 12 ] - Min[0i 3 ] unlike the previ- 
ous case because the contour of |ct| = tt does not reach 
at a upper-left corner of shaded region (this means that 
the point A does not correspond to Min[#i 3 ] ~ 7.0°). 
In this case, the point A where 6*13 ~ 8.0° determines 
Max[6» 23 ] as Max[<9 23 ] ~ 90° - Min[0 12 ] - 8.0° ~ 50.7°. 
It is also noticed from the FIG. [1] (b) that the value of 
Miners] ~ 7.3° is not determined at the point A but 
at C. Therefore, the value of Min[0i 3 ] ~ 7.3°, which 
leads to (9 23 ~ 90° - Max[6»i 2 ] - 7.3° ~ 45.2°, does 
not determine Max[# 2 3] as explained above. Regard- 
ing with Max(Min) [cr] , it corresponds to contours giv- 
ing Max(Min)[# 2 3] in the cases of m 3 = \J\ Am^ and 
VIAm^l+Aml!. 

We have also analyzed a case of 7713 = 0.2 eV. Re- 
sults are given in TAB. fl] In the case, values of |ct| 
within 102° < \a\ < 105° can give the solutions. Al- 
lowed region of mixing angles are the same as ones in 
the to 3 = yIAto|^J case, i.e. contours of solutions 
can rea ch at th e both points A and B like the case of 
m 3 = -\/|Ato|jJ . For the Max(Min) [a] , it is determined 
by Min(Max)[# 2 3] in contrast with the above two cases. 

We comment on running effects of CP-phases. It has 
been seen that the one of Majorana phases a is im- 
portant for the stability of self-complementarity. One 
may worry about running effect of a on low energy so- 
lutions for the stability, i.e. whether such effects spoil 
the solutions at a high energy scale or not. We have 
approximated the running effects of S, p, and a from 



the seesaw scale on solutions by using leading-log esti- 
mation in the RGEs. These running effects (A5, Ap, Act) 
are (AS, Ap, Act) ~ (0(0.001°), 0(0.001°), 0(0.001°)), 
(Q(l°), O(0.1°), O(0 1°)), an d (Q(l°), Q(l°), Q (l°)) for 
the cases of m 3 = ^\Am 31 \, ^/|A?«§ 1 | + Amfi, and 0.2 
eV, respectively. Therefore, most regions of a are stable 
against such small running effects up to the seesaw scale 
because of Max[tr] — Mui[ct] > Act in all cases. Further, 
small (but non-vanishing) running effects of S and p do 
not affect on the stability of self-complementarity The 
values in the table do not include such effects of devia- 
tions. Then one must remember that our analyses are 
only for the Majorana neutrino case. 

Next, let us consider case of \a\ = tt. It is worth em- 
phasizing that all running effects on CP-phases are van- 
ishing in the |ct| = tt case. The solution can be real- 
ized in narrow range as y/\ Am^ | + Amfi — 5.11 x 10~ 2 
cV< to 3 < 5.15 x 10~ 2 eV. The one of solutions of 
|ct| = 7r is shown in FIG. Q] (b). We find that points 
A and C constrain mixing angles to 31.3° < #i 2 < 37.5°, 
7-3° < 6»i3 < 8°, and Max(Min) [0 23 ] ~ 50.7°(45.2°) 
in the case of vanishing running effects of CP-phases 

(H=7T). 

We have also analyzed IH cases with to 2 = 
a/| Am 2 ^ I + Am 2 21 and 5.73 x 10" 2 eV. The value of 
ni2 = 5.73 x 10~ 2 eV is maximal one of m 2 giving so- 
lutions. These results are also shown in TAB. fl] In the 
IH case of m 2 = y^Am 2 ^! + Am 21 , there exist solu- 
tions in 78.8° < |ct| < 80.4°, and mixing angles are con- 
strained to 31.3° < X 2 < 37.5°, 7.3° < 13 < 11.1°, 
and Max(Min)[6> 23 ] ~ 51.4°(41.4°). For the m 2 ~ 
5.73 x 10~ 2 eV case, some solutions appear in narrow 
range as 87.6° < |ct| < 88.8°, 31.3° < 0x2 < 31.5°, 
10.7° < 6»i3 < 11.1°, and Max(Min)[6» 23 ] ~ 48.0°(47.4°). 

Finally, we evaluate effective mass term of neutrino- less 
double beta decay (0^/3/3), (m) ee = | SLi(^PMNs)e; m il; 
in our parameter space. It is written down as 

/ \ I 2 2 2ip 1 2 2 2ia , 2 I fo\ 

{m) ee = \mxc 12 c 13 e p + m 2 s 12 c 13 e +m 3 s 13 |, (8) 



4 



in our notation. The phenomenon of 0^/3/3 can distin- 
guish whether neutrinos are Dirac or Majorana particles. 
The results at benchmarks given in Tab. U are also pre- 
sented in the table. 

The magnitude of (m) ee strongly depends on the scale 
of mi or ?7i2 rather than mixing angles and CP-phases 
in the cases. In the NH with 7773 = y^Am|7l ancl 
y| Athjj I + Am| 1; dominant contribution comes from 
the second term of RHS of (|8|) because of the small sf 3 
and vanishing 777 1. In the cases of NH with 7773 = 0.2 eV 
and IH, the first term of (j8]) becomes dominant because 
of mass degeneracy for the NH and inverted mass 
spectrum for the IH, respectively. We predict 2.36 (14.9) 
meV< (m)ee < 97.4 (26.1) meV for the NH (IH) within 
the parameter space to make the self-complementarity 
stable. The Heidelberg-Moscow experiment [l(J for 
Qvj3f3 is giving the most severe bound on (m) ee , which 
is (m) ee < 210 meV. The CUORE experiment [Tl| 
is expected to reach (m) ee = (24 — 93) meV in the 
future. Therefore, a part of our predicting region may 
be checked in future experiments 

We have investigated stability of self-complementarity 
relation in the PMNS sector against quantum corrections 
by considering RGEs in the MSSM. We have found so- 
lutions stabilizing the self-complementarity both for the 
NH and IH cases by using low energy data of neutrino 
oscillation experiments. It has been seen that the one 
of Majorana phase a plays an important role to give the 
solutions. The self-complementarity relation can be sat- 
isfied up to an arbitrary high energy scale if neutrino 
parameters are correlatively within 31.3° < 6*12 < 37.5°, 
7-0° < 6>i 3 < 10.9°, and Max(Min)[6»2 3 ] ^ 51.7°(41.6°) 
with 126.0° < 1 0- 1 < 132.0° and 102° < a < 105° 
for the largest neutrino mass 7713 = Am^ | and 0.2 



cV of NH with 8 = p = at a low energy, res pec- 
tively. For the NH with 777,3 = \/l^ m 3il + ^ m 2ii we 
can obtain solutions in regions of 31.3° < #12 < 37.5°, 
7-3° < 0i3 < 10.9°, and Max(Min)[6>2 3 ] ^ 50.7°(41.6°) 
with 124.9° < \a\ < it, 5 = p = 0. For the IH case 

with the largest neutrino mass 7772 = -^/l^ 777 -!! I + ^ m 2i 
and 6 = p = at a low energy, the self-complementarity 
can be stable in regions of 31.3° < 6 12 < 37.5°, 7.3° < 
6»i3 < 11.1°, and Max(Min)[6»23] ~ 51.4°(41.4°). There 
exist solutions within Anijj < 7773 < 0.2 eV and 
•v/IAm^l + Ato^ < m 2 < 0.0573 cV for the NH and 
IH, respectively. In particular, the case of |er| = it and 
5 = p = 0, where the RGE running effects on CP-phases 
are vanishing, predicts 31.3° < 6 r2 < 37.5°, 7.3° < 
6»i3 < 8.0°, and Max(Min)[(9 23 ] ~ 50.7° (45.2°) for a nor- 
mal hierarchical neutrino mass spectrum (mi, 7772, 777,3) — 
(8.73 x 10~ 3 , 1.23 x 10" 2 , 5.11 x 10~ 2 ) eV. This kind of 
solution can be obtained in narrow range as 5.11 x 10 -2 
cV< 777,3 < 5.15 x 10~ 2 cV. 

Regarding with the 0^/3/3, the effective neutrino mass 
can be predicted as 2.36 (14.9) meV< (777} ee < 97.4 (26.1) 
meV for the NH (IH) within the parameter space to 
make the self-complementarity stable. The above solu- 
tions might be checked by combining more precise ex- 
perimental data of neutrino oscillations, determination 
of absolute neutrino mass scale, bounds on the sum of 
neutrino masses, and the discovery of 0^/3/3. 
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